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^^ ■ We study the projection of polytopes along k orthogonal vectors for various input 

and output forms. We show that if k is part of the input and we are interested 

O ' in output-sensitive algorithms, then in most forms the problem is equivalent to 

■^ ' enumerating vertices of polytopes, except in two where it is NP-hard. In two other 

forms the problem is trivial. We also review the complexity of computing projections 

p>^ I when the projection directions are picked at random. For full-dimensional polytopes 

containing origin in the interior, projection is an operation dual to intersecting the 

rj . polytope with a suitable hyperplane and so the results in this paper can be dualized 

r) \ by interchanging vertices with facets and projection with intersection. We would like 

^ ' to remark that even though most of the results in this paper do not appear to have 

O . been published before, they follow from straighforward reductions to other known 

results. The purpose of this paper is to serve as a reference to these results about 

the computational complexity of projection of polytopes onto affine subspaces. 

o 
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A polytope in M is a closed convex body that can be represented as either the convex hull of a 



f^ ' finite number of points or as the intersection of a finite number of halfspaces. We will call the 



former V-representation and the latter 7l!-representation. Accordingly a polytope will be called 
a V-polytope, an 7^-polytope or an 7-^V-polytope depending on whether the polytope is given 
by V, Ti. or both representations. For any polytope each of these representations is unique if 
^ ■ no redundancies are allowed and any of these representations completely determines the others. 

3 , We refer the reader to O E] for a thorough treatment of the subject. 

Suppose we are given an "H-polytope in M" and we want to compute the facets of its projection 
onto M"~^. One can use Fourier- Motzkin elimination (or any other algorithm of choice [7J) and 
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compute the facets easily. The projection can have as many as ^facets if P has m facets and as 
few as n facets. Now suppose we want to further project the shadow of P onto M"~^. Of course 
one can remove redundancies from the projection obtained in TZ"'^^ and apply Fourier-Motzkin 
elimination again but a careful reader might already have noticed a problem. If one wants to 
compute the projection of P onto a d-dimensional subspace, then applying Fourier-Motzkin 
incrementally might result in an intermediate polytope with very large (exponential) number 
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of facets even if we finally get a polytope with relatively small number of facets. Is there a way 
around this problem? 

This question is rather ill-defined at this point. To be precise, we want to know if there 
is an algorithm that computes the desired description of the projection and takes total time 
polynomial in the input and output. We need to include the size of the output in the picture 
because otherwise the problem clearly has no polynomial algorithm, since the number of facets 
of the projection of a polytope can be exponential in the size of the input. An algorithm 
that runs in time polynomial in the size of the input and the output is called output-sensitive 
polynomial algorithm. So the previous question can be rephrased as: Given an 7Y-polytope P, 
is there an output-sensitive polynomial algorithm that enumerates all facets of the projection 
of P onto some subspace? 

In this paper we answer this question in the negative. We also study several other natural 
variants of this problem arising from the simple fact that a polytope has two representations. 
One can ask the same question but in a different flavor by requiring the input and output to be 
in various different representations. We prove that in most forms this problem is equivalent to 
the problem of enumerating vertices of a polytope. In two forms, including the one described in 
previous paragraph, the problem is NP-hard and in two other forms the problem has a trivial 
solution. We also review the complexity of this projection problem if the projection directions 
are picked at random and prove that in many cases one can enumerate the desired form of the 
output in polynomial time. We would like to note that even though the algorithm for computing 
random projections was conceived independently by the author, it appears to be almost similar 
to the one presented in [8]. We nevertheless include it for completeness of discussion. 

The problem of enumerating the vertices of a polytope given by its facets has been studied 
for a long time by a number of researchers. Still the complexity status of Vertex Enumeration 
problem (VE), for general dimension and for polytopes that are neither simple nor simplicial, is 
unknown. It is neither known to be in P nor is it known to be NP-complete. The dual problem 
of computing "H-representation from V-representation is known as Convex Hull problem (CH). 
These two problems are equivalent modulo solving a Linear Program. Thus, for rational input 
these two problems are polynomial time equivalent and a polynomial output-sensitive algorithm 
for one can be used to solve the other in output-sensitive polynomial time. For more details 
about the problems with various Vertex Enumeration methods, we refer the reader to [1]. 

Our results about the hardness, and the equivalence of vertex enumeration and computing 
projection (in most forms) imply that in all forms where the projection can not be computed 
by a trivial algorithm, finding an output-sensitive polynomial algorithm will be a challenging 
task. Equivalently, an output-sensitive polynomial algorithm for vertex enumeration will have 
significant impact in many fields outside algorithmic polytope theory, like Control Theory, Con- 
staint Logic Programming Languages, Constraint Query Languages etc, where one frequently 
needs to solve the projection problem ([7]). 

In order to be able to talk about the equivalence of Vertex Enumeration and projection, we 
will define a complexity class based on Vertex Enumeration. Keeping in line with other no- 
tions of completeness, we call an enumeration problem <1> VE-complete if any output-sensitive 
polynomial algorithm for VE can be used to solve $ in output-sensitive polynomial time and 
vice-versa. Similarly, we call a problem VE-easy if it can be solved in output-sensitive polyno- 
mial time using an oracle for VE and we call a problem VE-hard if an oracle for this problem 
can be used to solve VE in output-sensitive polynomial time. 

The results in this paper are summarized in Table [U and Table [2l Table [U summarizes the 
complexity of computing projection along arbitrary directions while Table [2] summarizes the 
complexity of computing projection along randomly picked projection directions. 
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Table 1: Complexity of computing projection of a polytope onto an arbitrary subspace, for 
various input and output representations. 
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Table 2: Complexity of computing projection of a polytope onto a randomly picked subspace, 
for various input and output representations. 

Since for bounded polytopes containing the origin in relative interior, projection is an op- 
eration dual to intersection with a hyperplane passing through origin, all our results can be 
dualized by interchanging vertices and facets and replacing projection with intersection. 

The rest of the paper is organized as follows. In the next section we briefly review some related 
work about computing projection of a polytope. Our result section is divided into two parts. 
In Subsection 13.11 we present the results about the complexity of computing the projection of 
a polytope along arbitrary directions and in Subsection 13.21 we describe the complexity results 
for the case when the (orthogonal) projection directions are picked at random. 



2 Related Work 

Perhaps the best known algorithm for computing the facets of the projection of an 7i-polytope is 
the Fourier- Motzkin elimination discovered by Fourier in 1824 and then rediscovered by Motzkin 
in 1936. This method is analogous to the method of Gaussian elimination for equations and 
works by eliminating one variable at a time. Since eliminating one variable from a system of 
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m inequalities can result in ^ facets, the algorithm can have a terrible running time in bad 

cases where the intermediate polytopes have very large (exponential) number of facets but the 
final output has only a small number of facets. 

Many improvements have been done over the original algorithm (See [7] for a survey) but there 
is no algorithm that has an output-sensitive polynomial running time. The natural question 
then is whether one can find a shortcut around the intermediate projection steps in the Fourier- 
Motzkin elimination and obtain an output-sensitive polynomial algorithm. As we will see in 
Section [3l the answer is no if P 7^ NP. Thus the lack of any output-sensitive algorithm, 
for computing the facets of the projection of an "H-polytope, is somewhat natural because the 
problem turns out be NP-hard. 

Farka's Lemma provides a way of generating a cone whose extreme rays correspond to the 
facets of the projection of an 7i-polytope (j^), but unfortunately it does not yield a bijective 
mapping and many extreme rays of the resulting cone may correspond to redundant inequalities 
in the projection. Balas [3] found a way to get rid of these redundancies and provided a way 
to construct, given an T^-polytope P, another polyhedral cone W in polynomial time whose 



projection W yields a one-to-one correspondence between the extreme rays of W and the 
facets of the projection of P. Also, W has polynomially many facets compared to P. 

3 Results 

We will denote the projection of a n-dimensional polytope P onto a given d-dimensional subspace 
as TTdiP)- We will mostly omit the subscript and simply refer to the projection as tt{P) and the 
projection subspace will depend on the context. 

Formally, we are interested in the following problem: Given a polytope P S M" in 7Y-, V- 
or "HV-representation and a set F of fe orthogonal projection directions defining the projection 
space, we want to compute the non-redundant TL-, V- or TYV-representation of vr(P). 

3.1 Projection onto arbitrary subspaces 

Depending on the input and output form, we have nine variants of the problem. It is obvious that 
if the vertices are part of the input and one wants to compute the vertices of the projection, then 
each vertex can be projected trivially and the vertices that become redundant in the projection 
can be identified by solving one Linear Program per vertex. Hence, we have the following: 

Lemma 1. Given a polytope P C M" in V- or TiV -representation and an arbitrary projection 
subspace, non-redundant V -representation of it{P) can be computed in polynomial time. 

Also, it is easy to see that every polytope can be represented as the projection of a suitable 
simple^j- Moreover, given a polytope P by its vertices one can compute in polynomial time 
the vertices of this simplex A and the projection subspace such that P is the projection of A. 
Since it is trivial to compute the facets of a simplex given its vertices. Vertex Enumeration can 
be solved in output-sensitive polynomial time using any algorithm that computes the TC- or 
WV-representation of projection from V- or TYV-representation of a polytope. 

Clearly, one can also use any polynomial algorithm for Vertex Enumeration to compute the 
Tl- or TYV-representation of the projection of any polytope given in V- or TYV-representation in 
polynomial time. Hence, we have the following easy lemma: 

Lemma 2. Computing the TC- or TCV -representation of the projection of a polytope given in V- 
or TiV -representation is VE-complete. 

In what follows, we assume the input polytope P is of the form 

-P(x, y) = {x, y\Ax + By < 1} 

and we want to compute the projection 

7r(P) = {x e R'^\{x,y) G P for some y}, 

where A e Q™^'^,i? g Q"^^^^y £ M^. We also assume P to be full-dimensional and to contain 
the origin in its relative interior. For rational polytopes this assumption is justified because one 
can always find a point in the interior of the polytope via Linear Programming and move the 
origin to this point. We will sometimes omit details like A G Q'^^'^,B G Q^^'',x G M*^,?/ G M'' 
where it can be inferred from the context. 



^Assuming that the vertices are numbered through m ~ 1, simply append e^ to the i-th vertex, where eo is 
the zero vector and e^ is the i-th unit vector in R™^^. 



We are now left with the three cases where the input polytope is given by 7Y-representation. 
As we wiU see now, computing either the facets or the vertices of the projection in this case 
in hard while computing both facets and vertices of the projection is equivalent to Vertex 
Enumeration. Consider the following decision version of the problem: 
Input: Polytopes P = {x, y\Ax + By < 1} and Q = {x\A'x < 1} 
Output: YES if Q = 7r(P), NO otherwise. 

We will now prove that this decision problem is NP-complete thus proving the NP-hardness 
of the enumeration problem. 

Theorem 1. Given a polytope P = {x, y\Ax + By < 1} and Q = {x\A'x < 1} it is NP-complete 
to decide if Q j^ '^(P)- 

Proof. It is easy to see that deciding whether a given set of hyperplanes completely define the 
projection of a higher dimensional polytope is in NP. So the only thing remaining is to show 
that it is NP-hard as well. 

It is known ([5]) that Given an 7-^-polytope P = {x\Ax < 1} and a V-polytope Q = CHiV), 
it is NP-complete to decide whether P '^ Q. Clearly, P C Q if and only ii P Q = P. Now, 
P nQ has the following 7^-representation 

Ax < 1 

X — y^ Xv ■ V = 

v&V 

Xv > o,yv£v 

The variables A ensure that we consider only those points in P that can be represented as a 
convex combination of vertices of Q. One can further, in polynomial time, get a full dimensional 
representation of P fl Q by eliminating the d — 1 equations. The resulting polytope is a full- 
dimensional polytope in r!''^!"^. 

Since we are interested in only the (vector) variable x, the projection of this polytope along 
the axes corresponding to the variables A^, gives us the facets of PCiQ in the subspace of variables 
x. It follows that P f] Q = P if and only if the projection of the above defined polytope onto 
the subspace of variables x has the 7i-representation same as that of P i.e. Ax < 1. 

Thus any polynomial algorithm for deciding whether a given set of hyperplanes completely 
define the projection of some high dimensional 7^-polytope, can be used to decide whether an 
7^-polytope is contained in a V-polytope or not, which is a NP-complete problem. D 

Balas([3j) has shown that for a given ?-^-polytope (P) and a set of projection axes, one can 
compute the facets of another pointed polyhedral cone W and a set of projection axes such 
that the facets of vr(P) are in one-to-one correspondence with the extreme rays 7r(VF)G. The 
number of facets of W is polynomial in the number of facets of P. It is not difficult to modify 
the construction in [3j so that W is bounded i.e. a polytope and the vertices in the projection 
of W are in one-to-one correspondence with the facets of the projection of P. For completeness 
we state the result of Balas and describe the modification here. 



^Note that the projections of P and W are defined in different spaces and should not be confused as the same 
projection map despite the abuse of notation here. 



Lemma 3. Given an TC polytope P and a set of projection directions, there exists a polyhedral 
cone W and another set of directions such that the facets of tt{P) are in one-to-one correspon- 
dence with the extreme rays ofTT{W). Furthermore, W has polynomially many facets compared 
to P and the facets of W can he computed in polynomial time. 

We will use the notion of polar duality to prove that the cone W obtained from the construc- 
tion of Balas can be turned into a bounded polytope. For a polyhedral cone W in M" with 
facet inequalities Ax < and extreme rays the row vectors of V, where A and V are matrices 
with each row a vector in M". The polar dual W* has the roles of the extreme rays and facets 
reversed. In particular, the facet inequalities of W* are Vx < and the extreme rays are the 
row vectors of A. We again refer the reader to [H |^ for more details of the properties of this 
duality. 

Lemma 4. Given a pointed polyhedral cone W G M" and a set of projection directions T one 
can construct, in polynomial time, a polytope Q £ MJ^ such that the extreme rays of 7r(VF) are 
in one-to-one correspondence with the vertices of t:{Q), where both the projections are onto the 
same subspace. 

Proof. Clearly, none of the projection directions lie in the relative interior of W, otherwise the 
projection spans the whole subspace. Let W* be the polar dual of W. For any vector a in 
the relative interior of W* the hyperplane a ■ x = touches W only at the origin and hence 
W n {x\a ■ X < 1} is a bounded polytope. It is actually a pyramid with origin as the apex. 

Now consider the projection Tr{W), which is a pointed cone with origin as apex. This cone is 
a full-dimensional cone in the subspace containing it and we can consider its polar dual in that 
subspace. Let tt*{W) be the polar dual of tt{W). For any vector a in the relative interior of 
7r*{W), tt(W) n {x\a • x < 1} is a bounded polytope. Moreover, for such an a , 7j • q = for 
all 7j € r. Since tt{W) is a pointed cone such an a exists. 

Since it*{W) can be obtained as the intersection of W* with (^7 gf i^* ' 2; = 0}, a vector 
a in the relative interior of 7r*{W) can be computed in polynomial time if one knows either 
the vertices or facets of W. Also, a lies in the relative interior of W* as well. Define Q = 
Wri{a -x < 1}. Given the extreme rays (facets respectively) of W and the projection directions 
r one can compute the vertices (facets respectively) of Q in polynomial time. 

Since a is orthogonal to each of the projection directions, the vertices and facets of Tr{Q) are 
in one-to-one correspondence with the extreme rays and the facets of tt{W). D 

Theorem [1] together with Lemma [3] and Lemma H] gives the following: 

Theorem 2. Given a polytope P = {x,y\Ax -\- By < 1} and Q = CH{V) it is NP-complete to 
decide if Q ^ '^{P)- 

Now we consider the last variant of the projection problem where we are given an 7^-polytope 
and we want to compute the T^V-representation of the projection. As it turns out, although 
computing either the vertices or facets of the projection is NP-hard, computing both vertices 
and facets is VE-complete. 

Before we prove this, we would like to remark that the notion of output-sensitiveness can have 
various meanings. An output-sensitive polynomial algorithm for an enumeration problem (like 
VE) could enumerate vertices such that a new vertex is reported within incremental polynomial 
delay i.e. each new reporting takes time polynomial in the input and the output produced so 
far. It is equally conceivable that the algorithm takes total time polynomial in the input and 
output but there is no guarantee that successive reportings take only incremental polynomial 



delay. We will assume that if we have an output-sensitive algorithm of the latter kind, then we 
actually know the complexity of its running time. Under this assumption the two notions are 
same for VE. 

To see why this is true, consider the following. Given an 7i-polytope P and a V-polytope Q, determin- 
ing whether P = Q is polynomial time equivalent to VE (See [T]). Also, solving this problem gives an 
algorithm for VE that is not only output-sensitive polynomial but also has a polynomial delay guaran- 
tee. If we have an enumeration algorithm that has no guarantee of polynomial delay between successive 
outputs, but for which we know the running time, then we can use this procedure to create a polynomial 
algorithm for deciding the equivalence of 7i- and V-polytopes: Simply compute the time needed by the 
algorithm to enumerate all vertices of P assuming P = Q and run the enumeration algorithm for the 
time required to output |F| -I- 1 vertices. If the procedure stops then we can compare the list of vertices 
of P with that of Q in polynomial time. If, on the other hand, the procedure doesn't finish within the 
given time then P must have more vertices than Q and hence, P ^ Q. 

So for proving VE-completeness in the next theorem, when we assume the existence of an 
output-sensitive polynomial algorithm for VE, we also assume that this algorithm has a guar- 
antee of polynomial delay between successive outputs. Although we will work with the Convex 
Hull problem which is the dual version of VE, with a slight abuse of language we will refer to 
this dual problem as VE as well. 

Theorem 3. Given a polytope P = {x^y\Ax + By < 1} it is VE-complete to compute the facets 
and vertices of n{P). 

Proof. Since every polytope P G M" given by m vertices can be converted to a {m — 1)- 
dimensional simplex A such that P is a projection of A it is clear that computing TiV- 
representation of the projection of an T^-polytope is VE-hard. 

To prove that this problem is also VE-easy, we give an algorithm that uses a routine for VE 
to enumerate the facets and vertices of 7r(P). The algorithm proceeds as follows: At any point 
we have a list of vertices V of vr(P) and we want to verify that V indeed contains all vertices of 
7r(P). If the list is not complete, we want to find another vertex of vr(P) that is not already in 
V . To do this, we start enumerating facets of CHiy) and we verify that each generated facet 
is indeed a facet of vr(P). This is easy to check because of the following: 

Suppose h = {x\a-x = 1} be a hyperplane in the projection space. We say that h intersects P 

k times 

properly if the intersection P fl {(o, 0, ■ ■ ■ , 0) • (x, y) = 1} has some point in the relative interior 
of P. We will call such an intersection a proper intersection. 

We claim that the defining hyperplane of every facet / of CH{V), that is not a facet of vr(P), 
intersects P properly. To see this, pick a point xi in the relative interior of /. Such a point 
exists because CHiV) C vr(P) if some facet / of CHiV) is not a facet of vr(P). This point also 
lies in the relative interior of 7r(P). Also, there is a point (xi, yi) that lies in the relative interior 

k times 

of P that projects to xi. Clearly the hyperplane {(a, 0, • • • , 0) ■ (x, y) = 1} contains (xi, yi) and 
hence the hyperplane defining / intersects P properly. 

It follows that, if V does not contain all vertices of 7r(P) then there exists a facet / = {x|a-x = 
1} of CH[V) intersecting P properly. So if the enumeration procedure for facets of CH{V) 
stops and none of the facets intersect P properly then V contains all the vertices of vr(P). If 
some intermediate facet {a • x = 1} of CHiV) does intersect P properly then maximizing the 

k times 

objective function (a, 0, • • • , 0) • (x, y) over P produces a vertex of P that also gives a vertex v 
of ■k{P) upon projection. Moreover this vertex is not in the list V . Thus, if V is not a complete 



vertex description of tt{P) we can find another vertex of it{P) in polynomial time. This gives an 
output-sensitive polynomial algorithm for enumerating all facets and vertices of vr(P). Hence, 
computing all vertices and facets of the projection vr(P) of an 7^-polytope P is VE-easy as 
weh. D 

Given that it is widely believed that P ^ NP and VE has been studied quite closely by a 
number of researchers, one can infer that computing the projection is going to be a challenging 
problem for arbitrary projection directions. Also, one should note that if the projection is 
known to be a simple of simplicial polytope then computing both facets and vertices of the 
projection can be done in output-sensitive polynomial time because VE for simple or simplicial 
polytopes can be done in output-sensitive polynomial time ([2^,1]). But if the projection yields 
a degenerate polytope then we do not have such an algorithm. 

3.2 Projection onto random subspaces 

As stated before, the result that the facets of a random projection of a polytope can be computed 
in output sensitive polynomial time was independently established earlier by Jones et. al. (See 
[5). We nevertheless discuss random projections for completeness of the discussion and to 
highlight how the computational complexity of the projection problem changes when one is 
interested only in random projection directions. 

Now we describe a polynomial time algorithm for a more general case where the projection 
directions are, in some sense, non-degenerate. To make this notion precise, note that if P is 
the input polytope then every face of projection vr(P) is the shadow of some proper face of P. 
Call the maximal dimensional face /' of P a pre-image of the face / of vr(P) if / is obtained by 
projecting all vertices defining /' and taking their convex hull. In general, the dimensions of / 
and /' are not the same. This can happen if some projection directions lie in the affine hull of 
/'. We call a set of projection directions non- degenerate with respect to P if no directions lie 
in the affine hull of any face of P. 

Fact 1: If the projection directions are picked randomly, then they are non-degenerate with 
respect to any polytope P with probability 1. 

Also, for non-degenerate projection directions a face / of vr(P) and its pre-image /' have the 
same affine dimension. This is easy to see because a projection reduces the dimension of some 
face /' of P, that does not disappear in the projection, if and only if the projection direction 
lies in the affine hull of /' which is not possible for non-degenerate projection directions. Thus, 

Fact 2: For non-degenerate projection directions and a face / of 7r(P), if /' is the pre-image 
of / then dim{f) = dim(f'). 

Now, given a polytope P in ?Y-representation and a set of non-degenerate projection directions 
r we want to compute the facets of the projection '/r(P). Again, we assume that the facets of 
P are presented as inequalities of the form Ax < 1, where ^ is a matrix of size m x (d + k) 
and the projection has dimension d. Since, we will need to solve Linear Programs we also 
assume that the polytope is rational i.e. the entries in A are rational numbers. We will omit T 
from the discussion below and assume that the projection directions are aligned along a subset 
of co-ordinate axes. If not, we can apply a suitable affine transform to P depending on the 
orthogonal projection directions. Thus the reader should bear in mind that the polytope in 
what follows is a result of an affine transformation determined by a set of randomly picked 
orthogonal projection directions. 



Our algorithm for enumerating the facets of ir{P) proceeds as follows: Given a partial list of 
facets of 7r(P), for each facet / we identify its pre-image /' in P. For each of these faces of P 
we identify their (d — 2)-dimensional faces and among all such (d — 2)-faces of /' some give rise 
to ridges in tt{P). We identify which faces form the pre-image of some ridge of 7t{P) and from 
the corresponding ridge, we identify the two facets defining this ridge, thus finding a new facet 
of vr(P) if the current list of facets is not complete. 

Lemma 5. Given an TC-polytope P and a facet f of its projection tt{P), one can find the facets 
of P defining the pre-image of f in polynomial time. 

Proof. Let {x € M |a ■ x < l,a G Q } be the halfspace defining the facet / of vr(P). Clearly, 

k times 

the hyperplane h in M + with normal a' = (a, 0, • • • , 0) defines the supporting hyperplane 
{x € W^~^^\a' ■ X = 1}. Also, P Ci h is a face of P and is exactly the pre-image of /. A facet F 
of P contains this face iS P nhD F has the same dimensions as P Oh. Thus, one can find all 
the facets of P containing the pre-image of / in time polynomial in the size of P. D 

The next lemma follows immediately from the non-degeneracy of the projection directions, 
so we mention it without the proof (See Fact 2). 

Lemma 6. Given P and a facet f of its projection ir{P), if g is another facet of tt{P) sharing 
a ridge with f then the pre-images f and g' share a face in P. Furthermore, 

dim{f' n g) = dim{f Ci g) = dim{f') - 1 = dim{g') -1 = d-2 

Since the facets of P are known, we can identify all (d — 2)-faces of /'. The number of these 
faces is at most m for each pre-image /' and since /' is itself a polytope of dimension d— 1, we 
can compute the non-redundant inequalities defining the facets {d — 2-dimensional faces) of /o 
At this point, what remains is to identify these ridges and the facets defining these ridges. The 
following two lemmas achieve this. 

Lemma 7. Let P = {{x,y)\A ■ {x,y) < 1} be a polytope in R'^ x R'' , where A G Qrnx{d+k)^ 
X eR'^,y e R''. Also, let f be a {d-2)-face of P defined as f = {{x,y)\X ■{x,y) = l,{x,y) £ P}, 
where A C A. Then, f defines a ridge in the projection Q{x) if and only if 

k times 
pd 



there exists a € M such that (a,0, • • • ,0) S conv{A ), where each row of A is treated as 
a point in R'^'^^ . And, 



The feasible region of all such a is a line segment. 



It is not difficult to see that this lemma is just a rephrasing of the basic properties of sup- 
porting hyperplanes of a polytope. In other words, any hyperplane whose normal is a convex 
combinations of the normals of facets defining the face /, is a supporting hyperplane of P and 
vice- versa. Furthermore, if the normal lies in the subspace where the projection tt{P) lives, then 
it is also a supporting hyperplane of vr(P). Also, the normals of all hyperplanes that support 
a polytope at some ridge, when treated as points, form a 1-dimensional polytope i.e. a line 
segment. This formulation allows us to check in polynomial time whether a (d — 2)-face of P 
forms a pre-image of some ridge of Q{x). 



^Removing redundancies can be achieved via Linear Programming. 



Lemma 8. The end points of the feasible region of a in lemnia^ are the normals of the facets 
of Tr{P) defining the ridge corresponding to face f. 

As noted before, the normals of the hyperplanes supporting a polytope at a ridge r form a 
line segment when viewed as points. The end points of the segment represent the normals of 
the two facets defining the ridge r. This lemma ensures that given a pre- image of some ridge 
of tt{P), one can compute the normals of the two facets of 7r(P) defining the ridge r by solving 
a polynomial number of linear programs each of size polynomial in the size of input. 

Putting everything together we get the following theorem: 

Theorem 4. Given a polytope P defined by facets, and a set of non- degenerate orthogonal 
projection directions T one can enumerate all facets of it{P) in output- sensitive polynomial 
time. 

Since randomly picked projection directions are non-degenerate, one can also enumerate the 
facets of the projection of an 7^-polytope for such directions. Note, that this also gives an output- 
sensitive polynomial algorithm for the case when the input is an ?^V-polytope irrespective of 
the output form. Also, if the vertices of P are given then some tests like those in Lemma [7] 
and [H] become easier. We leave the proof of this to the reader since they do not affect our main 
argument about an output-sensitive polynomial algorithm. 

Corollary 1. Given an TCV-polytope P and a set of projection directions T that are non- 
degenerate with respect to P there is an algorithm that can enumerate the vertices and/or facets 
of 7t{P) in output-sensitive polynomial time. 

It is easy to see that computing the vertices of the projection of an 7-^-polytope along non- 
degenerate directions has VE as a special casqj. Hence, enumerating vertices of the projection 
of an 7^-polytope along non-degenerate projection directions remains VE-hard even though it is 
not clear if it remains NP-hard. Similarly one can argue that the complexity status of computing 
the projection of a V-polytope along non-degenerate projection directions remains the same as 
that of computing the projection along arbitrary directions. 
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